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The transport properties of a double quantum-dot device with one of the dots coupled to perfect 
conductors are analyzed using the numerical renormalization group technique and slave-boson mean- 
field theory. The coupling between the dots strongly influences the transport through the system 
leading to a non-monotonic dependence of the conductance as a function of the temperature and 
the magnetic field. For small inter-dot coupling and parameters such that both dots are in the 
Kondo regime, there is a two-stage screening of the dot's magnetic moments that is reflected in the 
conductance. In an intermediate temperature regime Kondo correlations develop on one of the dots 
and the conductance is enhanced. At low temperatures the Kondo effect takes place on the second 
dot leading to a singlet ground state in which the conductance is strongly suppressed. 
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Since the predictiorki*^ of the occurrence of the Kondo 
effect in a single quantum dot (QD) device and its subse- 
quent experimental observation^, several single and dou- 
ble quantum dot devices (DQDD) have been studied both 
theoretically'*''^ and experimentally.^'^ The interest in this 
systems stems form their potential applications to quan- 
tum and classical computing^iS and their usefulness as 
model systems for the study of the physics of strongly 
correlated electrons.® 

In its simplest form, the Kondo effect consists in the 
screening of a localized spin 1/2 magnetic moment anti- 
ferromagnetically coupled to a conduction-electron band 
and appeared originally in the context of magnetic im- 
purities in a metallic host. In a QD device Coulomb 
repulsion and spatial confinement result in charge quan- 
tization and its transport properties are dominated by 
the Coulomb blockade. When the charge on the QD is 
close to an odd integer the Kondo efect takes place and 
it results in perfect transmission through the system at 
sufficiently low temperatures in a 'wide range of gate volt- 
ages. 

When t"wo dots are coupled, a richer and more compli- 
cated behavior may arise^iSiiS including the possibility of 
observing quantum phase transitionsiii 

In this paper we present a detailed study of a DQDD 
in 'which only one of the dots is coupled to the leads (see 
FigureCJ. We discuss the case in 'which the QDs are very 
small and have a single relevant energy level at energies 
Eq and Eb and large charging energies Ua and [4, -where 
the indices a and b denote the t^wo dots. The position of 
energy levels of the dots and the strength of the coupling 
bet-ween them tab can be tuned by applying gate voltages 
in semiconductor devicesi^ 

The Hamiltonian of the DQDD is 



H — Hf) + He + Hn- 



D-E 



(1) 



'where 



e=a,b 



tab [di^dba + dl^daa^ 



(2) 




FIG. 1: Schematic representation of the device. A QD 
molecule with only one of the dots coupled to the leads. 



describes the a — b molecule and 

He = ea(k) 4^„Cko.a [a L,R) , (3) 

k,cr,a 

is the Hamiltonian of t'wo non-interacting source and 
drain leads. The coupling bet'ween the molecule and the 
leads is described by the last term in the Hamiltonian, 



Hd-e = 



dl 



aa kcra 



„t 



(4) 



For symmetric leads, the conductance through the sys- 
tem is-T 



G(T) 



doj ( pM, (5) 



'where /o^^(w) is the local electronic density of states on 
dot a. Here, A — 2npo{V^) 'where po is the electronic 
density of states of the electrodes at the Fermi level, the 
brackets denote the average over the Fermi surface, and 
/(w) is the Fermi function. 

We solved the above model using the numerical renor- 
malization group technique as 'well as an approximate 
analytical method that reproduces the main features of 
the numerical solution. Some exact results, that foUo'w 
from the application of Luttinger's theorem to this sys- 
tem, are also given. 

Our main findings are the following. The coupling tab 
is a relevant perturbation and the lo"w-energy physics of 
the DQDD is fundamentally different from that of a sin- 
gle QD device. The zero-temperature conductance of the 
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system vanishes when the total charge of the device is an 
even integer and it is perfect when the charge is an odd 
integer. This is an exact resuh. 

When both dots are nearly half-filled there are two 
different regimes for the temperature-dependence of the 
conductance depending on whether tab is small or large 
in a sense that will be made precise below. For small tab, 
a two-stage analog of the Kondo effect takes place with 
decreasing temperature. In the first stage, the magnetic 
moment of dot a is quenched by the electrons of the leads 
below some temperature. In the second stage, that oc- 
curs at a much lower temperature, the magnetic moment 
of dot b is in turn quenched by electrons that lie within 
the narrow Kondo resonance around the Fermi level. The 
two-stage Kondo effect results in a non-monotonic behav- 
ior of the conductance as a function of both gate voltage 
and magnetic field. In the regime of large tab, the mag- 
netic moments of the dots form a tightly bound singlet 
weakly coupled to the leads and the conductance varies 
monotonically at low temperature. 

When the charge on the device is varied by application 
of a gate voltage Vg, the zero-temperature conductance 
remains small in an extended interval of values of Vg. 
Conversely, in general, nearly perfect conductance is only 
observed in a very restricted range of values of the gate 
voltage. 

A surprisingly rich variety of possibilities exists for the 
dependence of the conductance upon gate voltage and 
temperature depending upon the various parameters of 
the model. 

The rest of the paper is organized as follows. In Sec- 
tion we derive exact results for the dependence of the 
zero-temperature conductance of DQDD upon its charge 
and magnetization. The solution of the model within 
slave-boson mean- field theory is presented in Section ITTl 
Section IIIII contains the description of our numerical re- 
sults. Finally, we summarize our conclusions in Section 

El 



I. EXACT RESULTS 

In the configuration depicted in Fig.^the ground state 
of a DQDD is expected to be a Fermi liquid. In this 
case, many exact results for the zero-temperature con- 
ductance of the device can be obtained using a simple 
generalization of Luttinger's theorem that we derive in 
the following. 

Since electron-electron interactions do not involve the 
electronic states of the leads, the latter may be integrated 
out from the start. All the properties of the DQDD can 
then be expressed in terms of the reduced 2x2 retarded 
Green-function matrix 

G-iH-Go-Hc.)-SH, (6) 

where Go(i^) is the non-interacting Green function, 

Go ^ 



UJ ~ Ea+iA tab 

tab — Cfc 



(7) 



and is the self-energy matrix. 

The total number of electrons with spin a in the dots 
at r = is given by 



duj Im [TrG(w)] 



(8) 



an expression that can be rewritten in the form 
=-ilmf° ^TrlnG-Vw) 

using the equality 



(9) 



TrG(a;) 



TrlnG"i(w) -KTr 



G(.)-S(.) 



(10) 

that can be easily checked using Eqs. jnj and Q. 

The second integral on the right-hand side of Eq. © 
vanishes order by order in perturbation theory in Ua^a^ 
which leads to 



1 



(11) 



where 



ip{uj) = Im In I [tj - ea - Saa(t^) + » A] 

eb - Sfcb(t^)] - [tab - Safc(w)]^ I • (12) 



Introducing the renormalized parameters ia = fa + 
Saa(O), €6 = £6 + Shb(O), and tab = tab - Sah(O) wc find: 



1 



<7 = 2 Imln [iah ~ ilb - iAh] 



(13) 



For a Fermi-liquid ground state the self-energy matrix 
at the Fermi energy is purely real and so are the renor- 
malized parameters. A straightforward calculation then 
yields: 



1 



Uc ^ — 



S , la and £(, > , ia h > ilf, 

27r + (5 , ea and h < , iah> ili, , (14) 

TT + (5 , otherwise 



where 



arctan 



^ah - ilb 



In all cases we have the relationship 



hA 



ab 



(15) 



(16) 



Expressions H14l) - (|16|l are identical to those that one 
would obtain for the non-interacting problem defined by 
Eq. ^ with the renormalized real parameters replacing 
the bare ones. 
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Similarly, the a-dot spectral density for spin a at the 
Fermi energy is given by 



pL(o) 



1 



\2 ' 



which, using Eq. 1)161) . can be cast in the form 
Paa(O) = ^sin^ (t™^) . 



(17) 



(18) 



Equations (0) and H18(l lead to the following expression 
for the conductance: 



5 -|; = ^E 



(19) 



where Go = 2e^//i is the quantum of conductance. 

In the absence of a magnetic field n1 — nj" and 
g = sin^ {ua +«{,)] where ne = ^^n^. The zero- 
temperature conductance thus vanishes when the total 
number of electrons in the DQDD is even. This occurs 
in particular when eg = —Ui/2 for both dots. At this 
point the Hamiltonian is electron-hole symmetric and 
Ua + Ub = 2. Applying a gate voltage Vg that shifts 
the levels of the dots, the charge of the device varies con- 
tinuously between Ua + Ub — (for Vg large and negative) 
and Ua + Ub = 4 (for Vg large and positive). The con- 
ductance g{Vg) vanishes at — >• ±oo and Vg = and 
has two maxima where it reaches the value g = I at the 
values of the gate voltage for which Ha + rib = 1 or 3. 

In the presence of a magnetic field B the dots polarize 
and a magnetization m — ^(n.| — ni) appears. In the 
regime Ua + Ub = 2 Eq. H19|l becomes 



g — sin^(7rTO). 



(20) 



In the low field limit g vanishes quadratically with B, 

g « TT^m^ « TT^B\^ = {B/B*f (21) 

were x is the spin susceptibility and B* = l/(7rx) defines 
a characteristic magnetic field. For very large B the dots 
become fully polarized, m — > 1, and the conductance also 
vanishes. At some intermediate field for which m ~ 1/2 
the conductance has a maximum at which .g = 1. 

II. SLAVE BOSON MEAN FIELD THEORY 



Here, Sf with £ = a,b are spin operators associated to 
the dots, and Sq = ^ J2s s' '4)s'^s,s'CQg, is the electron spin 
density on the orbital coupled to dot a. To leading order 
the coupling constant are Jab = '^t^bl i^a + Ub) and J = 
SA/nUa. 

Following standard method3i^*i^*i& we represent the 
spin operators as = ^J2ss'fes'^s,s'fis' in terms 
of two pseudo-fermions constrained by the conditions 
J2s flsfis — 1- The biquadratic interactions between the 
pseudo-fermions generated by the spin-spin interactions 
are decoupled introducing two Bosc fields, Ba (conju- 
gate to the amplitude fls'^os ^^^^ ^ab (conjugate to 
the amplitude flsfbs) ^^"^ ^^"^ constraints on occu- 
pations are implemented through Lagrange multipliers 
fields . The free energy expressed in terms of the Bose 
fields has a saddle point at which the latter condense, 

{Ba) = {fLcoa)^ {Bab) = {fLfaa)^ and the Lagrange 
multipliers vanish. The solution obtained by retaining 
only the saddle-point contribution to the free-energy is 
exact if the spin symmetry is extended from SU(2) to 
SU(Af) and the limit N oo \s taken. In the physical 
N — 2 case the SBMFT approach provides a simple and 
yet qualitatively accurate description of the low-energy 
properties of the system. 

At the saddle point the effective Hamiltonian is 

^^cff = JabY.{fUb. + h.c.) + jY.'^fLco. + h.c.) + HE, 

(7 a 

(23) 

where the renormalized couplings are determined self- 
consistently by the equations 



where 



Jab 
1 



D 



duj f{uj) Im Gabi^^) 



(24) 



-D 



duj f{uj)Regaa{co) , (25) 



-D 



ijj +JA -Jab 
— Jab ^ 



(26) 



and D is the half-width of the conduction band. 

The system of equations (|24|l - (|26|) can be easily solved 
numerically. However, we shall focus below on some rel- 
evant limiting cases for which analytical expressions can 
be obtained. 



In this Section we present a mean-field slave-boson 
treatmentiii^ (SBMFT) of the problem. For simplic- 
ity we restrict ourselves to the electron-hole symmetric 
case in the most interesting regime in which the Coulomb 
repulsion on the dots Uafi ^ A. In this case the occupa- 
tion of each dot is ~ 1 and the low-energy excitations 
of the DQDD are spin fluctuations. These are described 
by the effective Hamiltonian 

Hk = Jab^a ■ Sb + JSq • So + He- (22) 



A. Decoupled dots 

For Jab = Jab — the model reduces to the well 
known Kondo Hamiltonian. The results for this model 
in the SBMFT approximation are well known. In the 
weak-coupling limit Jpo ^ 1 the solution of Eq. H25|) 
at T = is Ao = De^yP"'^ . At finite temperatures 
solutions with non-zero Ag only exist below a tempera- 
ture Tc = ciDe~^/''°'^ ^ where ci is a constant of order 
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1. This approximation thus gives a spurious transition 
at which the spins and the conduction electrons decou- 
ple. In the exact solution of the Kondo model this sharp 
transition is replaced by a crossover at the Kondo tem- 
perature T% oc De-'^/P"-' . We thus identify with T%. 
The quasiparticle density of states is 



Ao 



(27) 



describing a resonance of width Aq at the Fermi level. 
The presence of this resonance leads to perfect conduc- 
tance, g — 1, at zero temperature. The dot's magnetic 

susceptibility is Xaa oc 1/ ^tt Aq j . 

B. Weakly coupled dots 

In the limit in which the coupling between the dots is 
weak, Jab <C Ao , the solutions of the equations at T = 
to leading order in Jab / A are 



Jab 



A 



Aexp 



ttA 

"iJab 



Ao 



1-2 In ^ 

V Ao / I Jab , 



(28) 



(29) 



While A « Ao remains essentially unchanged in this 
regime, the effective coupling between the dots is strongly 
suppressed. With increasing temperature Jab further de- 
creases and vanishes at a temperature 



To = C2Aoexp(- 



ttAo 



J' 



(X 



Jab ' Ao 

where C2 ~ 1 depends weakly on the parameters. As 
above, this transition temperature should be interpreted 
as a crossover temperature below which the magnetic mo- 
ment of the dot b is screened. 

Note that To -C T^ which suggests a two-stage screen- 
ing process. First, the magnetic moment of dot a is 
screened by the electrons of the leads below T^. Then, 
at the much lower temperature To, the magnetic moment 
of dot b is screened by the heavy quasiparticles of the lo- 
cal Fermi liquid that forms on dot a for T ^ T^. The 
form of To supports this physical picture as it corresponds 
precisely to the expression of the Kondo temperature of 
a magnetic moment screened by electrons of a band of 
width ^ Aq and density of states l/(7rAo). 

The quasiparticle densities of states of the dots are 



(30) 



Paa(w) 



Pbbi(^) 



1 



-Ji 



;A 



JahA 



'"[^'-Jlb) +(-a)" 



(31) 



(32) 



respectively. 

The first of these equations determines the conduc- 
tance through the device. In the temperature range 
To ^ T ^ T^, where Jab is irrelevant, we recover the 
usual Kondo resonance 



1 



A 



This leads to a large conductance g ~ 1- 
For T < To and < J^b, 



(33) 



Paa(t^) 



1 



^Ao ^2 



(34) 



Therefore, paa vanishes at the Fermi energy where it 
has a dip of width cx Tq. This dip is the analog of the 
Kondo hole that appears in the conduction-electron den- 
sity of states in the usual Kondo problem. An immedi- 
ate consequence of the presence of this hole is that the 
zero-temperature conductance of the DQDD vanishes at 
T = 0. 

In the same limit, the magnetic susceptibility of the 
DQDD is dominated by the contribution of dot b. We 
find 



X « Xbb 



/2 1 

ttAq To 



(35) 



Therefore, the magnetic scale B* defined in Eq. ()21|l is 
also determined by Tq. 



C. Strongly coupled dots 

For Jab > Aq there are no solutions of Eqs. 1)24(1 and 
((25|l with non- vanishing A. The dots remain decoupled 
from the leads and Jab = Jab/ '2- The ground state is 
a singlet and there is a energy gap Eg — 2 Jab = Jab 
in the electronic spectrum. Note that the singlet-triplet 
gap that results from the SBMFT treatment is the exact 
result for a pair of isolated dots. 

The decoupling of the dots from the leads is an arti- 
fact of SBMFT. Presumably, the fluctuation corrections 
neglected in this approach generate a weak coupling be- 
tween the spin singlet and the leads but we have not 
checked this point explicitly. 



III. NUMERICAL RESULTS 

In this section we present the numerical solution of the 
DQDD model obtained using the numerical renormaliza- 
tion group method (NRG)i2iiSiiS modified to improve the 
accuracy in the computation of the spectral density:^ We 
choose the half-bandwidth of the conduction band D as 
the unit of energy and consider for simplicity the case of 
identical dots, £a = £b = £ and Ua = Ub = U. 
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We shall first discuss the electron-hole symmetric case 
in order to make contact with the SBMFT results of the 
previous Section. 



A. Kondo screening 

Figure |21 shows the NRG results for the temperature 
dependence of square of the total magnetic moment of the 
DQDD in the electron-hole symmetric case for several 
values of tab- The other parameters are U = 0.5, e = 
-0.25, and A = 0.035. 

The magnetic moment ^ is defined through the re- 
lationship = Tx where x is the contribution of the 
quantum dots to the total magnetic susceptibility of the 
systemjii For a free S=l/2 spin, fi^ = 1/4. 

The uppermost curve in Fig. El represents the case of 
decoupled dots, tab = 0. At high temperatures (but low 
compared to U) the spins of both dots are free. Then 
/i^ — > 1/2. At some temperature scale there is a drop in 
/i^ due to the Kondo screening of dot a. We identify this 
temperature with the Kondo temperature (a more 
precise determination of will be given in the following 
Section). Since the spin on dot b remains free down to 
T = 0,fj,'^ = 1/4 for T < TO . 

For any finite tab the total magnetic moment in the 
ground state vanishes and fi^ — > as T ^ 0. All the 
curves but the lowermost correspond to values of tab such 
that Jab < T^. It can be seen that, within this range of 
parameters, screening takes place in two stages as pre- 
dicted by SBMFT. With decreasing temperature there is 
a first drop in /i^ at the scale T'^ followed by a plateau. 
A second drop occurs at a much lower scale that we shall 
identify below with Tq. The temperature at which the 
first drop occurs is rather insensitive to the value of Jab 
in agreement with Eq. H29|l . Note that the magnetic mo- 
ment varies logarithmically with temperature at the two 
steps which confirms the two stages of screening are of 
the Kondo type as stated in the previous Section. 

When tab increases the steps become closer, which sig- 
nals a merging of the two scales with increasing Jab, in 
agreement with Eq. H30|l . 

We found that, in this low- Ja& regime and at low tem- 
perature, all the curves collapse onto a universal curve 
if the temperature is scaled by an appropriate factor Tab 
that depends on Jab- We expect on physical grounds 
Tab oc Tq. In order to check this hypothesis we calcu- 
lated Tab from the susceptibility data using the criterion 
fJ-^{Tab) = 1/8. The results are represented with solid 
circles in the inset to Fig. [21 It is seen that Tab varies 
exponentially with t~^ in agreement with Eq. H30|l . The 
solid line is a fit to the SBMFT expression for Tq. An 
arbitrary multiplicative constant in the definition of Tab 
was absorbed in C2 and the parameter ci is the propor- 
tionality factor between T'^ and Ag. As expected, the 
values of Ci and C2 determined by the fit are 0(1). 

The lowermost curve shown in Fig. (21 corresponds to 
the case Jab > Tk- seen that it is qualitative differ- 




FIG. 2: Magnetic moment of the DQDD as a function of 
the temperature for different values of the interdot coupling. 
From top to bottom, tab ~ 0, tab = 0.002 — 0.0046 in steps of 
0.0002, and tab = 0.02. Parameters are U = 0.5, e = -0.25 
and A = 0.035. Inset: binding energy of the a — b singlet Tab 
estimated using the NRG (see text). The solid line is a fit 
based on Eq. (QUI . 



ent from the others. In this regime the spins of the dots 
form a robust singlet weakly coupled to the leads. The 
susceptibility and the magnetic moment decrease expo- 
nentially in Jab/T. 



B. Spectral properties and conductance 

Figure O shows the spectral density on dot a at zero 
temperature in the symmetric case for U — 0.5 , A = 
0.035 and tab = 0.0022. This value corresponds to the 
third curve from the top in Fig. (21 

We observe three features, a very narrow central peak 
and two broad Coulomb peaks located at w = ±J7/2. 
We determine the Kondo scale T^ = 7 x 10^^ from the 
full width at half-maximum of the central feature. The 
insets to the figure are zooms of the central peak. The 
left inset shows the Kondo peak plotted as a function of 
the reduced frequency lj/T^. The right inset shows the 
Kondo hole in the density of states present over the tiny 
frequency interval \uj\ < Tq = 1.78 x 10~^°. This feature 
is obviously not visible on the scales of either the main 
plot or the left inset. 

The conductance of the DQDD is entirely determined 
by the spectral density of dot a [ci. Eq. jSJ]. The tem- 
perature dependence of g in zero magnetic field is repre- 
sented in the lower panel of Fig.^lfor tab — -ix 10""^ with 
the rest of the parameters fixed at the values given above. 
For temperatures T > T^ the device is in the Coulomb 
blockade regime and the conductance is low. With de- 
creasing temperature the Kondo correlations start to 
build up and the magnetic moment of dot a decreases 
as shown in the upper panel of the figure. For T < T^ 



6 



10 



s 



2 







10 
8 












10 

8 




1 






1 ' 










1 ' 1 ' 1 






3 6 

CL 4 












3 " 

Qr4 


- 








2 












2 













, 1 , 











1 , U , 1 






5 


-2.5 2.5 


5 


-4 


-2 2 


4 






— ^ 












r— ~ — 





FIG. 3: Spectral density of dot a at T = 0. Parameters are 
U = 0.5, A = 0.035 and tab = 0.0022. The broadened atomic 
levels of dot a are clearly seen a.t uj = ±U /2. The left inset 
shows the central Kondo peak of width T]^ ~ 7 x 10~*. The 
right inset shows the Kondo hole of width To ~ 1.78 x 10"^". 
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FIG. 4: Upper panel: Magnetic moment of the DQDD at 
zero field as a function of temperature for tab = 0.003. The 
other parameters are as in Fig|5| Lower panel: Conductance 
through the device at zero field as a function of the temper- 
ature. Three regimes can be clearly seen both in the conduc- 
tance and the magnetic moment. 

the spin on dot a is quenched resulting in resonant scat- 
tering through dot a and an enhanced conductance. For 
temperatures T < Tq ^ 7 x 10~^, ho"wever, the spin of 
the second dot becomes also screened and the Kondo hole 
appears in the density of states. In this regime the con- 
ductance decreases again. 

We no"w turn to a description of our results in the pres- 
ence of a magnetic field. Data arc displayed in Fig. [S] for 
the same parameters as in Fig. 0] The upper panel sho"ws 
the field-dependence of the square of ^-component of the 
total spin S = Sa + Sb of the DQDD at T = 0. The 
lo"wer panel sho"ws that of the conductance. 

For fields B ;» the spins of both dots are fully 




10" 10""* 10"" 10° 10" 



FIG. 5: Upper panel : magnetic moment of the DQDD as a 
function of magnetic field at T = 0. Lower panel: conduc- 
tance as a function of the magnetic field at low temperatures. 
Inset: Field-dependence of the T = conductance. The cou- 
pling between the dots tab = 0.003. The other parameters are 
as in Fig|21 

polarized and {SD — 1. The Kondo effect can not take 
place and the conductance is small. 

For B < the first stage of the Kondo effect takes 
place and the spin on dot a is quenched. At this point 
{SD ~ 1/4 indicating that the spin on dot h still remains 
polarized by the field. In this region of fields the con- 
ductance increases "with decreasing field, approaching its 
maximum value, g — 1. 

Upon further decrease of the field, Kondo screening of 
the second spin starts and, for B <To ^ 10^'^ x T^, it is 
in turn fully quenched. Then, the conductance decreases 
again and vanishes as i? — > 0. 

The inset to the figure sho"ws that the lo"w-field conduc- 
tance varies quadratically -with B "with a characteristic 
field B* ~ To in full agreement "with Eqs. (|2U and 

The effect of the temperature on the field-dependence 
of the conductance is also sho'wn in the figure. For T = 
10"'' T°. w To and fields B > T the results are very 
similar to those obtained at zero temperature. At lo"wer 
fields, ho"wever, the conductance saturates. This behavior 
arises because, at this temperature T >Tq, quenching of 
spin b is only partial and the Kondo hole in the density 
of states ceases to develop. For T > both stages of 
the Kondo effect are suppressed by thermal fluctuations 
and the conductance is featureless as a function of B. 



IV. CONDUCTANCE AS A FUNCTION OF 
THE GATE VOLTAGE 

In this Section we discuss the dependence of the con- 
ductance upon the gate voltage. For simplicity we re- 
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strict the analysis to the case of identical dots Sa = Sb = 
e, Ua = Ub = U. We have checked, however, that the 
results for other situations are qualitatively similar. 

We display in Fig. IHlthe conductance as function of e, 
calculated from Eq. 0, for the set of parameters U = 
0.25, tab = 0.025 and A = 0.125. In this case the Kondo 
temperature and the effective magnetic coupling between 
the dots are r^- ~ A and Jab ~ 0.01 < for e = -U/2. 

At T = we observe two peaks of perfect conductance, 
g = 1, and a wide region of low conductance between the 
two. The conductance vanishes at e = —U/2 as discussed 
in the previous Section. 

We have checked that the conductance at T = sat- 
isfies Luttinger's theorem by computing g from Eq. H19|) 
using the charge of the DQDD obtained from the NRG 
calculation and comparing it with the outcome of the di- 
rect calculation. The two results are the same within the 
numerical uncertainty. 

The peaks of perfect conductance g = 1 occur when 
the DQDD is in either of the charge states Ua + nb = 1 
or Ua + rib = 3. At large e the charge of the DQDD is 
very small and increases with decreasing e. It reaches 
n = 1 near e* such that the lowest energy level of the 
a — b "molecule" crosses the Fermi level. In the absence 
of interactions this is e* — tab- Interaction effects will 
reduce e*, an effect that is already present for weak in- 
teractions for which a simple Hartree-Fock aproximation 

leads to s* = ijf^,^ + (US/A)'^ — U/A, where S — nb — Ua is 
the difference between the charges on the two dots. For 
our paramenters the reduction is very large and we find 
e* ^ 0.001 < tab- 

For temperatures < T < T^, the conductance in- 
creases at the center of the valley due to the disappear- 
ance of the second stage of the Kondo efect. This leads 
to the three peak structure observed in the second and 
third panels from the bottom. If the temperature is fur- 
ther increased, the height of the conductance peaks de- 
creases and, for temperatures T ~ A, there is a single 
broad conductance peak. For our values of the parame- 
ters there are no Coulomb blockade peaks because their 
width A is comparable with their separation U. 

Figure shows the function g{e) at several tempera- 
tures for a different set of parameters, Ua — Ub — 0.5, 
tab = 0.05, and A = 0.063. In this case, ~ 0.01 and 
Jab ~ 0.02 > T^. This a case in which the physics at 
the symmetric point is dominated by the formation of a 
strongly bound singlet between the spins of the two dots. 
The behavior at zero temperature is qualitatively similar 
to that found in the previous case. The evolution of the 
conductance with temperature is however quite different 
on two accounts. First, the conductance at e = —U/2 
is always a minimum; second, two additional conduc- 
tance peaks appear upon rasing the temperature in the 
region < T < T^. This phenomenon can be under- 
stood by noticing that while both Jab and increase 
with £ in the valley away from e = —U/2, the increase of 
is more pronounced because, in this case, the depen- 
dence on £ appears in the exponential. As result, while 




FIG. 6: Conductance through the double dot system as a 



function of the gate voltage Vg 



-£ and 



different temperatures. Parameters are Ua — Ut — 0.25, tat = 
0.025 and A = 0.125. 




FIG. 7: Same as in Fig.|S|but for a different set of parameters: 
Ua^Ut^ 0.5, tab = 0.05 and A = 0.063. 



Jab > Tk ^'^^ e = there is an inversion of the inequality 
for sufhciently large e. When this inversion takes place 
the normal Kondo effect can occur and results in an en- 
hancement of the conductance around some intermediate 
point. Finally, at high temperatures the usual Coulomb 
blockade peaks are this time visible because, for this pa- 
rameter set, A ^ [/. 

We close this Section with an analysis of a case in 
which the hopping matrix element between the dots is 
large. Fig. |H1 shows results for tab = 0.25 = U/2 with 
all the other parameters as in Fig. \7\ The appropriate 
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smaller or larger than the Kondo temperature of a single 
dot. When Jab ^ T'^ the transport properties of the sys- 
tem reflect the existence of a two-stage Kondo effect. In 
the first one, the magnetic moment of the dot connected 
to the leads is screened at , the Kondo temperature of 
the isolated dot. In this regime the conductance increases 
with decreasing temperature. At temperature Tq <^ 
the magnetic moment of the second dot is also quenched, 
a process that leads to a conductance decreasing with 
temperature. A similar non-monotonic behavior is ob- 
served in the field-dependence of the conductance. In 
the opposite case. Jab 3> T^, the magnetic moments of 
the dots bind forming a strong singlet weakly coupled to 
the leads and the low temperature conductance decreases 
monotonically with T. 



FIG. 8: Same as in Fig. |7| but for tat = 0.25. The curves 
correspond to temperatures ranging from 9.2 x 10~^T^ to 
18.4 from bottom to top (in the central region). Each 
curve corresponds to a temperature that is twice as high as 
that of the previous one. 



starting point for a qualitative description of this case is 
an isolated a — b molecule whose bonding and antibond- 
ing orbitals have energies e± = e ± tab- Consider the 
case in which the energy of the bonding orbital crosses 
the Fermi energy. The probability of occupation of the 
antibonding orbital is now strongly reduced because of 
the large energy gap 2tab that separates its energy from 
that of the bonding state. States with non-zero occupa- 
tion of the antibonding state can then safely be projected 
out of the subspace of available states using a general- 
ized Schrieffer- Wolff transformation just as one projects 
out states with double occupancy of either of the dots. 
The resulting effective problem describes a single quan- 
tum dot with renormalized parameters. We thus expect 
to observe the ordinary single-dot Kondo effect around 
£ — tab and e = —tab — U. This is precisely what the 
figure shows. At low T, there are two peaks of width 
^ U in the conductance in the range of values of e for 
which the occupation number of the DDQD is odd. In 
the peak region wc observe Coulomb blockade and the 
Kondo effect depending on the temperature. ^ 



V. CONCLUSIONS 

In summary, we studied the spectral, magnetic and 
transport properties of a double quantum-dot device in 
which one of the dots is coupled to perfect conducting 
leads. 

The zero-temperature conductance of the system van- 
ishes when the total charge of the device is an even integer 
and it is perfect when it is an odd integer. 

When both dots are individually in the Kondo regime 
two different physical situations are possible depending 
on whether the magnetic coupling between the dots is 



Application of a gate voltage leads to a rich variety of 
behaviors for the dependence of the conductance upon 
gate voltage and temperature. A general feature is that 
the conductance is small in a very wide range of values 
of Vg at all temperatures. At low temperatures, regions 
of nearly perfect conductance do exist, but they are re- 
stricted to narrow intervals of values of the gate voltage. 
Phenomena similar to some of those reported here have 
been found in a study of a multilevel quantum dotiSi 

In a recent paper— (see also Ref. l23j) the same model 
was studied using exact diagonalization methods in small 
clusters. These authors investigate the zero-temperature 
conductance of the device as a function of the positions 
of the dot's energy levels ea,b- In the regime of small tab 
they report perfect conductance for Sa = £b = —U/2. 
These numerical results are inconsistent with ours and, 
most importantly, with the implications of Luttinger's 
theorem. The exact diagonalization method apparently 
fails to capture the second stage of the Kondo effect re- 
sponsible for the supression of the conductance at zero 
temperature. 

The reason for this failure is easy to understand. In 
order to correctly describe the single dot Kondo effect 
in a finite cluster, the level spacing in the leads must be 
smaller than the Kondo temperature— iSSiSLSLSS For a 
QD coupled to a linear chain the number of sites required 
is Ng 3> 2D/Tk which becames numerically intractable 
beyond Ng ^ 20. Large Kondo temperatures must then 
be considered. The Kondo temperature of the second dot 
is much lower than Tk rending the problem intractable 
by the exact diagonalization method in the limit of small 
tab- However, these zero-temperature exact diagonaliza- 
tion results describe qualitatively the behavior for tem- 
peratures To < T <C or for large tab- 

Finally, we want to stress that our conclusions are valid 
in the case in which the effective magnetic coupling be- 
tween the dots Jab is antiferromagnetic and do not apply 
directly to the case of ferromagnetic coupling. 
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